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Abstract 

By constructing the general six-parameter bright two-soliton solution of 
the integrable coupled nonlinear Schrodinger equation (Manakov model) us- 
ing the Hirota method, we find that the solitons exhibit certain novel inelastic 
collision properties, which have not been observed in any other (1+1) dimen- 
sional soliton system so far. In particular, we identify the exciting possibility 
of switching solitons between modes by changing the phase. However, the 
standard elastic collision property of solitons is regained with specific choices 
of parameters. 

PACS numbers: 42.81. Dp, 42.65.Tg, 03.40.Kf 

1 



Recent developments in the field of optical solitons in fibers have already revealed 
the possibility of overcoming the limitations on the speed and distance of linear- 
wave transmission systems [0,0. However, the interaction between optical solitons 
influences directly the capacity and quality of communication [Q, ||. Since the 
nonlinear Schrodinger (NLS) equation is used as a mathematical model for such 
studies, considerable attention has been devoted to analyze the nature of collision 
between NLS solitons |], |J §, £§ . In more general physical situations || , coupled NLS 
equations are often used to describe the interaction among the modes in nonlinear 
optics, as for example, in the case of birefringent || and other two-mode |7j fibers. 
Here we consider the integrable coupled NLS equation of Manakov type M, 



(1) 



where q% and q<i are slowly varying envelopes of the two interacting optical modes, 
the variables x and t are the normalized distance and time, and \i is a positive 
parameter. 
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Exact soliton solutions have been derived |], [H| |TT| for the system (|l|) with dif- 



ferent procedures. In |J, using bright one-soliton solutions of the system ([!]), many 
physical phenomena such as birefringence property, soliton trapping and daughter 



wave ('shadow') formation are studied. Further, as noted in ||12|| , the Manakov model 
([!]) is important in describing the effects of averaged random birefringence on an or- 
thogonally polarized pulse in a real fiber. In addition, by considering the analytic 
solution of the system ([]]), conditions have been established for soliton switching and 
energy coupling among the two modes in a nonlinear fiber [D3|. Recently in [ID|, two 



of the present authors (R.R and M.L) have derived bright and dark multi-soliton 
solutions of the system ([!]). The interaction between bright and dark solitary waves 



is studied in 14 



In this paper, we focus our attention on the system (|1]) to derive a more general 
bright two-soliton solution than the ones that have been presented in the literature 
[ |T0|j by using the Hirota method [13], IB|. The asymptotic behavior of this solution 



is studied further in order to explain the collision properties of coupled solitons in 
the anomalous dispersion region. 

It is often stated in the literature || [17| that the two-soliton solution is very 
difficult to derive and too complicated to analyze. [Nevertheless Manakov || was 
able to make an asymptotic analysis of the solitons.] In this paper we present for 
the first time the completely general two-soliton solution to ([!]) in explicit form, 
which can be readily used asymptotic analysis, numerical computations, etc. Our 
explicit solution not only clarifies Manakov's observation on the asymptotic behavior 
but also helps to realize the exciting possibility of a novel type of inelastic collision 
allowing switching between components via phase change; this in contrast to the 
standard elastic collision usually observed in (l+l)-dimensional soliton systems. 
[However, the standard elastic collision property of the solution is recovered when 
restrictions are imposed on some of the free parameters.] The situation is reminiscent 
of the dromion solutions in (2+l)-dimensional systems such as the Davey-Stewartson 
equation |T9f where inelasticity has been observed [^, |2(| in the scattering 



process. The details of the present study are as follows. 

Recently, a special form of the bright two-soliton solution for Eq.(ffl) with five 



arbitrary complex parameters has been derived in |10j, using the Hirota method 



[ |15| , |T6fl . The first step in this method is to transform the system (|I|) into the Hirota 
bilinear form. For this purpose we use the transformations [[l(|, q\ = g/f and 
<?2 = h/f, with / real, to obtain the bilinear form of (jl|) as 

(iD x +D*)g-f = 0, (iD x + Dl) h ■ f = 0, D 2 J ■ f = 2/i (gg* + hh*) . (2) 

Here the operator D is defined by D™a ■ b = (d Xl — 9 X2 )a(xi)6(x 2 )| a;=Xl=X2 . The 
one-soliton solution to (§) is given by / = ae m , h = /3e m , g = 1 + e m+v ^ +R , that is, 



ae 



Pe<» r ^(|a| 2 +|/?| 2 



where 

r]j = kj (t + ikjx) + r]f . (4) 
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The arbitrary complex parameters in (|3|) are k\ and any two of the set a, (3, rji . 
This solution can also be written in the more conventional form 

cosh(r] 1R + <p) 



where 



\ A 1 = %, A 2 = ^ A = ^/z(M 2 +l/?l 2 )- 



r A , A , 

and we have introduced the subscripts i? and J for the real and imaginary parts of 
the quantity in question. (A positive root has been used to define ez R and hence 
is real.) From this form it is easy to identify the amplitude A and the phase 0. We 
also note that km and ku determine the amplitude and velocity of solitons, and in 
Manakov's notation || A4 corresponds to the unit polarization vector of the soliton. 

The two-soliton solution can be obtained by substituting into (0) the expansion 

9 = X9i + X 3 93, h = xh + X 3 h, and / = 1 + x 2 f2 + X 4 fa, where \ is the formal 
expansion parameter. The main problem is to choose the proper ansatz for g 1; h±. 



In [[TIJ a bright two-soliton solution was derived in this way assuming the input 
expression g% = e Vl + e m , hi = e £ (e m +e m ). The resulting two-soliton solution 
shows the standard elastic collision with a phase-shift. 

However, for any integrable equation it must be possible to combine any pair of 



one-soliton solutions into a two-soliton solution ||16|| , and therefore it should also be 
possible to start with gi and h\ given by 

9l = ai e^ + a 2 e m , h x = (3^ + f3 2 e m , (6) 

and this way generate a more general bright two-soliton solution with six arbitrary 
complex parameters k\, k 2 , ax, a 2 , j3\, and (3 2 (Note that the parameters T]j in 
(4) have been absorbed into aj and (3j). This is indeed possible. By following the 



equations (10-28) in [10), we obtain the most general expressions for the two optical 
modes q\ and q 2 as 



I _|_ e vi+Vi+Ri 4- e m+V2+ 5 o _|_ e v*+v2+$o _|_ e v2+V2+ R z -f e vi+v*+m+ , n2+ R 3 ' 



(7) 
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+ (3 2 e m + e Vl+v ^ +V2+5 ' 1 + e Vl+V2+ri 2 +s! 2 

where [note that Ri are real] 

+ fei + jfej;' k 2 + k* 2 ' {) 

e&2 = J^V^^W)^ 12 '^ 22 ^ (10) 

e " = ®^ (12) 

e * 3 = (fc 1+ ^(t+"^ 1+ ^l 2(K11K22 " K12K21) ' (13) 

and = ^ + Ag& . (M ) 

/b^ ~r n^- 

Does the introduction of the additional parameters make any qualitative change 
in the behavior of soliton? The answer is yes and we find the novel result that the 
above general solution (0) corresponds to an inelastic collision of two bright solitons, 
as long as a\ : a 2 ^ (3\ '■ /3 2 . In order to see this, we analyze the asymptotic form of 
the solution (0). 

One important advantage of Hirota's method is that the solution allows easy 
analysis of the asymptotic behavior by taking a given rjm — > ±oo in ([?[) and com- 
paring the result with the one-soliton solution (Q). The interpretation of the result 
in terms of the actual motion of the solitons depends on the signs of km and ku. 
In general r] iR = k iR (t — 2kux) so that soliton j is located in the vicinity of the line 
t = 2kjjx. Let us change to the frame comoving with soliton j (coordinatized by £) 
by putting x = (t — £)/(2fcjj). Then rjj R = kj R £, while for the other soliton m we 
get r] mR = k mR (l - ^f)t + k mR ^f£. Thus if we have, for example, k iR > and 
k\ik 2 i < 0, which corresponds to a head-on collision, we find that rjm — > ±oo corre- 
sponds to t — > ±oo. [Correspondence between the signs of x — > ±oo and r] iR — > ±oo 
is obtained in a similar way and leads to a slightly different dependence on the signs 



of the fe's.] 

In each limit rjm — > ±00 the resulting q% can be written as in (H) with different 
amplitudes A and phases 0. Let us denote by A\ the amplitude of the component 
qi of the soliton k as the other soliton goes to ±00, and similarly for the phase <fi. 
Furthermore let us define the phase shift <3> fc = <ft k+ — <f) k ~, and a "transition matrix" 
T k by A k+ = A k ~T k . [Since the magnitudes of A and T are the most interesting 
quantities we will not study their phases here.] We find the following results: 



Soliton 1: 



\A\-\ = \a 1 \/A 1 , \Af\ = \p 1 \/A 1 , (15) 
x |1 - A 2 a 2 /ai| ■ M _ |1 - Ag^/Aj / 1R x 

1 — A1A2 V 1 — A1A2 



ii?!, ^-(i^-i?!-^)- (17) 



Soliton 2: 



AM = |a 2 |/A 2 , \Af\ = |/3 2 |/A 2 , (18) 

1^2, |l-Aiai/a 2 | ,^2, |1 - Ai/3i//3 2 | 
Mil - ? =- M2 1 - ~7= ; l iy J 

y |1 — AiA 2 | y 1 1 — A1A2I 

= \B,2, $ 2 = i(i?3-i?l- J R 2 )(=$ 1 ). (20) 



where 



Aj = \l^{\a,i\ 2 + |/5j| 2 ), Ai = K2i//cn, A 2 = ki 2 /k 22 . (21) 

The expressions above are scale invariant in the way that, for example, allows us to 
take j3i = 1. If we also have ci\ = a 2 then one can easily verify that |T/| = 1, which 
implies perfect elastic scattering. We also note that the noninteracting stationary 



pulses discussed in [17] follow with the choice a.% — fl\ — 0, ku = kn = 0. 



One can verify that \A\~ | 2 + |A 2 ~| 2 = \A\ + \ 2 + |A 2 + | 2 = 1/fi so that the total 
intensity of each soliton is conserved. However, the distribution of this intensity 
among the two component fields can change during collision. It turns out that an 
inelastic effect can be obtained just by changing the relative phases of the parameters 



In figure 1 we have a picture of a head-on collision with hi — 1 + i, k 2 — 2 — i, 
Pi — 1, cei = 1 and a 2 = (39 + i80)/89. We still have \a 2 \ = 1, but its nontrivial 
phase (a 2 = e te , 9 xs 64°) is enough to cause quite dramatic non-elasticity. The 
initial time-profiles at both ends of the x-axis are evenly split between the two 
components, but the profile observed later at the large positive x end is almost 
completely in component 2. If we had chosen a 2 — 1 the scattering would have been 
elastic. Thus, by changing the relative phase of just one component of one soliton 
we are able to change the final state quite dramatically. 

Finally, we wish to connect explicitly our results on inelastic collision with that 
of the work of Manakov in [8]. In formulating the results that he obtained by an 
asymptotic analysis of the inverse scattering problem associated with (1), Manakov 
used a (complex unit) polarization vector obtained from the amplitudes of the soli- 
tons ((19) and (20) in [8]). Correspondingly, in our case the initial polarization 
vectors are 



a,- 



^Af/VkP + IAI 2 , (22) 

while the final vectors can be expressed as 

1 h-h / k 2 + fc* \ 

Cl = xhTwA c '-WT¥j c >- c >^)' (23) 

^ = i rrM*-rTF<* < 24 > 

X k 2 + kf \ k 2 + kf J 

where x = WTWlV \hTW lci '° 2 1 ) • (25) 

(The connection between Manakov's ( and our k is £i = (—ik 2 )* : ( 2 = —ik 2 .) 

Manakov pointed out that during soliton collision their velocities and amplitudes 
(intensities) do not change but the associated unit polarization vectors do change 
provided they are neither parallel nor orthogonal. Our observation, illustrated in 
Fig. 1, is that even if we keep the initial amplitudes unchanged and just change their 
relative phases we can change the amplitude distribution in the final polarization 
vector. The parametric choice associated with Fig.l is c\ = {\/y/2, l/y / 2) T , c 2 = 
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(-(39 + 80i)/(89 v / 2), l/v / 2) T ,Ci = 1 + C2 = -1 - 2i. Each component here 
has magnitude l/y/2. After collision the unit polarization vectors take the values 
£/ = (-0.0190 + 0.0593^, -0.734 + 0.677i), c 2 ' = (0.422 + 0.272i, 0.462-0. 731i), with 
magnitude distribution (0.0623, 0.998) and (0.502, 0.865), respectively. One can also 
note that the magnitude of the first component of t\ can be even made exactly zero 
if we can also change the parameters hi suitably. However, in practical applications 
it is probably only the relative phase that can be easily modified. 

To conclude, the general two-soliton solution (|7|) of the (l+l)-dimensional sys- 
tem ([I]) exhibits a novel type of inelastic collision, not seen in any other (1+1)- 
dimensional evolution equation. Naturally such a property will have important ram- 
ifications in optical fiber communication such as providing intensity pump sources, 
soliton switching and so on, which remain to be explored. 

It will also be of interest to investigate the ramifications of the above type 
of inelastic collision in non-integrable cases, for example when the nonlinear cross 
coupling coefficients are different from one. 

The work of R.R. and M.L. forms part of a Department of Science and Tech- 
nology research project. The work of J.H. is partially supported by the Academy of 
Finland, project 31445. 

Figure Captions 

Fig.l Intensity profiles |gi| (top) and \q 2 | (bottom) of the head-on collision solution 
(|7D with the parameter values k\ = 1 + i, k 2 = 2 — i, fa = {3 2 = oci = 1 and 
a 2 = (39 + i80)/89. 
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Inelastic Collision and Switching ... by Radhakrishnan et al 
Figure 1 




